Preliminaries:
library(tidyverse)
library(ggplot2) # for pretty plots
library(dplyr) # for %>%
library(infer) # for statistics
library(tidyr) # nice tables that look nice
library(car) # ANOVA
library(formatR) #to format R
library(MASS) # reported used for normalizing data for ANOVA
[Include a brief summary of the paper you are reanalyzing data from (e.g., the overall objective of the paper, the types of data collected and how sampling was done, what the main results were) and lay out what you intend to replicate.]
[Include a view of the first few lines of the dataset plus any exploratory data analysis - e.g., overview of descriptive statistics included in the paper, plots, etc.]
load in dataset:
f <- "https://raw.githubusercontent.com/slcornett/data-analysis-replication/main/data/Serrano_et_al_2020_MPMH.csv"
d <- read_csv(f, col_names = TRUE, show_col_types = FALSE) # show column names, hides dataframe message details
head(d)
## # A tibble: 6 × 19
## Nombre Captura Sex Sexo Track Calls5min Temp HR Larvas LHC Peso
## <chr> <chr> <chr> <chr> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 "Eugenia" 142 Fema… H 376L 8 NA NA NA 23.4 1.15
## 2 "Michelle" 158 Non-… M 162L 4 15.4 84.2 NA 22.4 NA
## 3 "Dr\xe1cul… 171 Fema… H 172R 5 NA NA NA 23.7 1.1
## 4 "Herlinda" 173 Non-… M 171L 5 NA NA NA 22.7 0.92
## 5 "Luisa" 175 Fema… H 179L 11 NA NA NA 23.2 1.06
## 6 "Femalesra… 176 Preg… MP 174L 3 NA NA 5 22.6 1.11
## # … with 8 more variables: SPL <dbl>, NC <dbl>, ND <dbl>, CD <dbl>,
## # `Agg Entropy` <dbl>, `Avg Entropy` <dbl>, DF <dbl>, `RMS Amp` <dbl>
DESCRIPTIVE VARIABLES: Nombre = name given to ID the frog <3; sexo = sex shorthand (MP = pregnant male, H = female, M = non-pregnant male); Captura = number captured Track = the call recording track
HR =
Larvas = number of babes in their mouth! PHYSICAL MEASUREMENTS: Peso = weight/mass (g); LHC = longitud hocico a cola (snout-vent length, mm); Temp = internal frog temperature (ºC)
ACOUSTIC VARIABLES: Calls5min = calls in 5 min interval (call repetition rate); SPL = Sound Pressure Level (dB); NC = number of notes/call; ND = note duration (ms); CD = call duration (s); Agg Entropy = aggregate entropy; Avg Entropy = average entropy; DF = dominant frequency of call (Hz)
RMS Amp = Root Mean Square Amplitude
[Be sure to thoroughly explain what replications you are doing and comment your code so that it is easy for a reader to understand. Include in this section relevant tables/figures/values from the original paper for comparison to what you accomplished with your replication. Note that I want you to do the bulk of any exposition using text and markdown syntax outside of code blocks. That is, your document should not just be one big code block with R style comments but rather a nicely formatted report with code separated from exposition, interpretation, and dicussion.]
Narrative section that overviews how successful were you at replicating the analyses and visualizations in the study.
1. What problems did you encounter?
2. Why might you have encountered those problems?
3. What details were lacking from the original study’s methods that might have hampered your ability to replicate the authors’ results?
[Include the citation for your paper, plus any other literature you might reference.] Serrano, J. M., Penna, M., Valenzuela-Sánchez, A., Mendez, M. A., & Azat, C. (2020). Monomorphic call structure and dimorphic vocal phenology in a sex-role reversed frog. Behavioral Ecology and Sociobiology, 74(10), 127. https://doi.org/10.1007/s00265-020-02903-3